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We study the SIR epidemic model with infections carried by k 
particles making independent random walks on a random regular 
graph. We give a edge-weighted graph reduction of the dynamics 
of the process that allows us to apply standard results of Erdos- 
Renyi random graphs on the particle set. In particular, we show how 
the parameters of the model produce two phase transitions: In the 
subcritical regime, 0(ln k) particles are infected. In the supercritical 
regime, for a constant C G (0, 1) determined by the parameters of the 
model, Ck get infected with probability C, and O(lnfe) get infected 
with probability (1 — C). Finally, there is a regime in which all k 
particles are infected. Furthermore, the edge weights give information 
about when a particle becomes infected. We demonstrate how this 
can be exploited to determine the completion time of the process by 
applying a result of Janson on randomly edge weighted graphs. 

1. Introduction. The spread of an infection throughout a population, 
often referred to loosely as an epidemic, has come to be modelled in var- 
ious ways in the literature, spurred by the richness of domains in which 
the abstract notion of a virus has gone beyond the traditional biological 
phenomenon. Electronic viruses over computer networks are not the only 
extension, others include rumour spreading [22] or broadcasting [3] and viral 
marketing [2]. Models may vary over domains, but the underlying principle 
is one of spread of some unit of information or state through interaction 
between individuals. 

In much of the literature on the spread of epidemics as well as the dis- 
semination of information, individuals reside at fixed vertices of a graph 
and the evolution of the state of an individual depends on the state of its 
neighbours in the graph. In particular if the graph is complete, mean-field 
(homogeneous mixing) models have been exploited to study the outcome 
of diffusion process [9]. More recently, there has been an increasing interest 
in understanding the impact of the network topology on the spread of epi- 
demics in networks with fixed nodes, see [13] for a review of such results. 
There has however been little analytical work to date on models where the 
possible interactions between the nodes are dynamic, i.e., evolves in time. 

We explore a particular instance of dynamic interaction by assuming that 
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individuals are mobile particles and can only infect each other if they are in 
sufficiently close proximity. The model is motivated both by certain kinds of 
biological epidemics, whose transmission may be dominated by sites at which 
individuals gather in close proximity (e.g., workplaces or public transport 
for a disease like SARS, cattle markets for foot-and-mouth disease, etc.) 
and by malware. Furthermore, it is relevant to studying the dissemination 
of information in opportunistic networks [5] where the information is trans- 
mitted between users who happen to be in each others range. As in the case 
of static networks [22] one may be interested in the time it takes for the 
rumour to be known to all users. 

In our model (elaborated upon below) there are k particles making inde- 
pendent, discrete-time, synchronous random walks on an n-vertex r-regular 
random graph G. Each particle is in one of three states: susceptible (S), 
infected (I), recovered (R). An infected particle can infect a susceptible par- 
ticle, which remains infected for a fixed infectious period before recovering 
permanently. This is known as the SIR epidemic model and is extensively 
studied in the literature. When ^ = oo (the SI model) particles never go 
from I to R. 

Two questions can be asked: (1) When ^ < oo, how many particles ever 
get infected? (2) When ^ = oo how long does the process take to complete? 
We address both of these questions by reducing the dynamics of the process 
to what we term an interaction graph. This turns out to have the structure 
of an Erdos-Renyi (E-R) random graph Qk,q on the set of particles, where 
the edge probability g is a function of the parameters of the model. Infected 
particles are connected components in Gk,q, and so well known results from 
the literature on E-R random graphs can be directly applied using our re- 
duction to answer question (1). In particular, we show how the parameters of 
the model produce two phase transitions: In the subcritical regime, 0(lnA:) 
particles are infected. In the supercritical regime, for a constant C deter- 
mined by the parameters of the model, Ck get infected with probability C, 
and O(lnfc) get infected with probability (1 — C). Finally, there is a regime 
in which all k particles are infected. Statements are with high probability 
(whp), that is, with probability tending to 1 as n — t- oo. 

Furthermore, the interaction graph reduction assigns weights on the edges 
that give information about when a particle becomes infected. This infor- 
mation can be used for addressing question (2). As an example, in the case 
of /9 = 1,^ = oo, we apply a result of Janson [16] on randomly-weighted 
graphs to deduce that completion time, the time at which the final infection 
takes place, converges in probability to Ink, where 6^ = jE^- 

Whilst the metaphor of an epidemic is a motivating and illustrative one, 
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this work is part of the more general scope of interacting particle systems 
(see, e.g., [1, ch. 14]). 

2. Model. Let r > 3 be fixed. Let Qr denote the set of r-regular 
graphs with vertex set V = {1,2, and the uniform measure. Let 

G = {Vg,Eg) be chosen uniformly at random (uar) from Qr- The results in 
this paper are always asymptotic in n = \ Vg\- The notation An ~ Bn means 
that lim„_!.oo An/Bn = 1- The notation O, o, 0, have their usual meanings 
with respect to n. We denote by f^^(l) a quantity that can be replaced by 
an arbitrarily large constant. 

At step let S{t),I{t),Tl{t) be the set of susceptible, infected and recov- 
ered particles respectively. These form a partition of the particle set V. Let 
/o = |X(0)|. 

When two particles x and y meet at some vertex v at time step t, an 
interaction takes place between them at that time step with probability 
p £ (0, 1], which is a fixed constant parameter of the model. We term such 
an event an xy interaction and call p the interaction probability. If one of 
the particles is infected and the other is susceptible, the infection takes place 
upon interaction. The infectious period is either cxd or a positive integer. 

Consider that the time step counter has just changed from t — ltot: every 
particle x makes an independent move in its random walk. Subsequently, the 
rules are as follows: 

1. If X, y are on the same vertex v, there is an xy interaction with prob- 
ability p ; if they are on different vertices they cannot interact ; each 
particle pair (non)interaction is independent of every other particle 
pair (non) interaction 

2. If X G S{t — 1) and if there is an xy interaction with some y € I{t — 1), 
then X G I{t), otherwise, x G S^t); in the former case, we say x was 
infected at time t 

3. If X was infected at time tx, then x G T^{t) if tx + S, < t, otherwise 
X G I{t) 

(note, we assume that 7^(0) = 0, and so a particle could only be in TZ{t) if 
it had been infected at some time previous to t). 

Observe two things from the above rules. Firstly, infections are not tran- 
sitive in a single time step. For example, suppose x G I{t — 1) meets 
y,z G S{t — 1) at vertex v at time step t. If x interacts with y but not 
z, then y does not pass on the infection to z at that time step, regardless of 
whether or not they interact. 

Secondly, the rules allow, in principle, for infections to be passed on even 
with infectious period = 1; there could be a chain of infections with x 
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infecting y, which in turn infects z, etc. 

A note on notation and terminology: When we refer to a 'period' [ii,t2]) 
we mean the set of time steps {ti, ti + l,ti + 2 . . . , ^2}- When we refer to 
"time t" we are referring to step t on the counter - this is a discrete time 
process. The first time step is t = 1, and t = is not a time step, but a 
useful convention to express the state of the system before the first time 
step. 

3. Assumptions. We first specify some assumptions of the model. 

If each particle is at distance at least cj(/c,n) = Q{\n\nn + ln k) from every 
other then we say the particles are in general position (g.p.). We assume the 
following: (i) G is typical, (ii) The number of particles k <n'' where e > 
is some sufficiently small constant, (iii) Particles start in general position, 
(iv) /o = |X(0)| > 1 and |7^(0)| = 0. 

A graph G is typical if it satisfies certain conditions. The definition will 
be given in following sections but for now it suffices to say that most graphs 
in Qr are typical, that is, a graph G picked uar from Qr will be typical whp. 

Assumption (iii) is also not unreasonable; it is straightforward to verify 
that if the positions of each of the k particles are chosen uar from the 
vertices of G, then whp they will be in g.p. with respect to each other if e 
is small enough. 

4. Results. We state the results for the case Iq = 1, that is, there is 
initially a single infected particle. It is relatively straightforward to generalise 
to Jo > 1, and our proofs cater to this. 

Let Mk be the total number of particles that ever get infected in the 
course of the process, and let be the completion time of the process; the 
time step at which the last infection takes place. Define 

(1) ^=^^^-^) 



-2 + p 

(observe that p <ip <V). 

Theorem 1. Suppose fc — t- 00 as n 00. Define 

(2) ^ = k(^- 

where Or = Then 

(i) If(f)<l then whp, = 0(lnA;) + 1. 
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(ii) If (j) ^ c where c > 1 is some constant, then with probability (l + o(l))C, 

Mfc = (1 + o{l))Ck, otherwise = 0{lnk) + 1. C is the unique 
solution in (0, 1) of the equation 1 — x = e~^^. 

(iii) If (j) > (1 + e)liik where e > 1 is some constant, then whp M^, = k 
(i.e., all the particles get infected). 

Theorem 1 is for finite ^, but observe that taking the convention that 
= for < 1 means that part (iii) is consistent with the SI model, for 
which all particles get infected almost surely. The theorem effectively gives 
conditions for transitions between different possible "regimes" of behaviour. 
The most interesting is the regime of part (ii), which is entered as (p transi- 
tions from <j) < 1 to (j) > 1. When /q = 1, roughly speaking, in this transition 
the number of infected particles goes from very few (0(lnA:)) whp, to a 
constant fraction Ck with probability C, or 0(lnA;) with probability 1 — C. 
Hence, it's "all or nothing", with a constant probability C of being "all". 

Concerning the completion times, we will demonstrate in section 12 that 
for ^ = oo, (that is, the SI model) how the edge weightings can be ex- 
ploited by a fairly straightforward application of a theorem of [16] to get the 
following. 

Theorem 2. 



where T^. is the completion time for k particles. 

5. Related Work. In this section, we briefly describe some of the rel- 
evant related work on diffusion processes like epidemic spreading and the 
dissemination of information in mobile environments. 

There has been a growing body of work in the interacting particle systems 
community analysing epidemic models with mobile particles. In [10] the 
authors provide a review of the results, techniques and conjectures when the 
graph is an infinite lattice. In [23], the authors explore by means of mean- 
field approximations the evolution of the number of infected individuals 
when individuals perform random motions on the plane. Recent papers by 
Peres et al [21], Pettarin et al [24] and Lam et al [19] analyse mobile networks 
modeled as multiple random walks; as Brownian motion on in [21], as 
walks on a 2-dimensional grid in [24] and as walks on a grids of dimension 3 
and above in [19]. In each case, there is a parameter r within which distance 
a pair of walks can communicate, producing a communication graph (which 
is disconnected below a certain threshold Tc). [21] studies various aspects of 



(3) 



nlnk/k 



Tk 
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the communication graph, such as how long it takes a particular point to 
become part of the giant component. [24] studies the broadcast time Tb of 
a piece of information and originating in one agent in relation to r. Setting 
r = means information is passed only when particles are coincident. In this 
case, Tb is equivalent to our completion time and [24] give, for k > 2, Tb = 
@{n/\/Tc) whp^. In [19] r = 1 and the results they have show a significant 
difference to the 2-dimensional case of [24]. We state their results for the 3- 

1 2 - 

dimensional case: There exists a constant c such that if ens log < k < n 
then Tb = e{ni/Vky, if k < cn3 log ^713 then Tb = @{n/k). Statements 
are whp. 

Of closer relevance to this work are [11] and [12]. Both of these papers 
study infections carried by random walks on graphs. In particular [11] anal- 
yses an SI model similar to ours; multiple random walks on a graph G that 
carry a virus, and the quantity of interest is the completion time. They give 
a general upper bound E[rfe] = 0{m* Ink) for any graph G, where m* is the 
expected meeting time of a pair of random walks maximised over all pairs of 
starting vertices. Special cases are analysed too, in particular, they give an 
upper bound of E[rfc] = In Inn) for random r-regular graphs. This is 
a factor Inn larger than the precise value of the process considered in this 
paper. 

Finally, we mention [3], which studies flooding on dynamic random net- 
works. A fixed set of n vertices is part of a dynamic random graph process 
where each edge is an independent two-state Markov chain, either existing 
or not existing. A single initial vertex initially holds a message and any 
vertex which receives this message broadcasts it on existing edges for the 
next k steps. Although flooding is a different process to multiple random 
walks, the authors develop a reduction to a weighted random graph with 
some similarity to the interaction graphs we present. It allows them to de- 
rive relations between the edge-Markov probabilities and state asymptotic 
bounds on the number of vertices that receive the message, as well as the 
broadcast (equivalently, completion) time. 

6. Outline of Paper. In the next section, we lay the foundation for 
our analysis by defining typical graphs and describing aspects of random 
walks on them. We then give a series of lemmas which serve as tools for our 
proofs. 

Following that, we deal first with the case p = 1 for both SI and SIR 
models. Section 8 treats the SI model, and defines the interaction graph. In 

^The tilde notation hides polylogarithmic factors, e.g., 0(f{n)) = 0(/(n) log"^ n) for 
some constant c. 
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it we calculate the probability of getting a particular outcome of the process, 
that is, a particular interaction graph. 

Section 9 builds on section 8 to take into account a finite infectious period. 
This will give us Theorem 1 for the special case p = 1- 

We then move on to the general case p < 1. To do so, we first demonstrate 
the general principles and techniques of our approach by looking at a system 
of only two particles (section 10). This allows us to convey the intuition 
behind the more involved calculations given in section 11 which deals with 
k particles. This section is a p < 1 generalisation of the SI model of section 
8. The case of p < 1,.^ < oo follows immediately since it extends section 11 
in precisely the same way that section 9 extends section 8. 

Finally, in section 12 we deal with the completion for the SI model in the 
general case /? < 1 to get Theorem 2 

7. Approach. In this section we detail key concepts and lemmas that 
we use to analyse the viral process. In section 7.1 we give a formal definition 
of typical graphs. Knowing the properties of typical graphs, we can make 
statements about how walks behave on them. In section 7.2 we describe 
how the long-term behaviour of a random walk - specifically its convergence 
to the stationary distribution - relates to the eigenvalues of its transition 
matrix. We also show how to map the k multiple walks on G on to a single 
walk on another graph, the product graph H. Much of the analysis of meeting 
times between walks is done through the framework of the product graph. A 
particular pair of particles meeting in G maps to the single walk on H being 
at a set of vertices of H. Such a set of vertices is contracted to a single vertex, 
resulting derived graph T upon which a single walk moves. Thus calculations 
of particle meeting times in G is done by calculating the hitting times of a 
single walk to a contracted vertex in T. Lemma 4 shows that for each of the 
graphs G, H, T, the walks are rapidly mixing, meaning that they converge to 
their respective stationary distributions quickly. This is a crucial component 
of the proofs and the behaviour of the processes. In section 7.3 we introduce 
a key lemma. Lemma 5, which allows us to make precise calculations of 
T-u), the probability that a walk starting on u visits v for the first time 
(after mixing) at time t. When this lemma is applied to contracted vertices 
in r, it gives us probabilities for meeting times between particles walking 
on G. Lemma 7 in section 7.4 gives formal justification for the use of vertex 
contraction to reason about visits to sets and section 7.5 investigates the 
probability that a particular vertex v was visited given that a set of vertices 
S with V £ S has been visited. 

Lemma 5, in conjunction with Lemmas 7, 8 and 9 culminate in Lemma 
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10 in section 7.6. This lemma gives probabilities for a particular pair of 
particles meeting at time t and no other pair meeting before hand. It is the 
main tool used in section 8. 

7.1. Typical Graphs. We say an r-regular graph G is typical if it satisfies 
the properties P1-P4 hsted below. We first give some definitions. 

Let 

(4) Li = [ei log,, nj , 

where ei > is a sufficiently small constant. 

A vertex v is treelike if there is no cycle in the subgraph G[v, Li] induced 
by the set of vertices within (graph) distance Li of v. 

A cycle G is small if |C| < Li. 

PI. G is connected and not bipartite. 

P2. The second eigenvalue of the adjacency matrix of G is at most 2^/r^^^+ 

e, where e > is an arbitrarily small constant. 
P3. There are at most n^^^ vertices on small cycles. 

P4. No pair of cycles Gi,G2 with |Ci|,|C2| < lOOLi are within distance 
lOOLi of each other. 

Note that P3 implies that at most n^"^ vertices of a typical r-regular graph 
are not treelike, where 

(5) n'^ = O (r^in^^i) = O (n^^^) . 

Lemma 3 ([8]). Let Q[. C Q^. he the set of typical r-regular graphs. Then 

\Q'r\^\Qr\- 

PI implies that a random walk will converge to a stationary distribution 
TT on the vertex set. Because the graph is regular, vr will be the uniform 
distribution. P2 implies that a random walk will converge quickly to the 
stationary distribution; it will be rapidly mixing. In fact, for all the graphs 
we consider, O(fclnn) steps is sufficient for our results. A typical graph also 
has most of its vertices treelike. 

7.2. Convergence to stationarity and product graph formulation. Let G 
be a connected graph with n vertices and m edges. For random walk Wu 
starting at a vertex u of G, let Wu{t) be the vertex reached at step t. Let 
P = P{G) be the matrix of transition probabilities of the walk, and let 
Pu{v) = Pr{yVu{t) = v). If the random walk Wu on G is ergodic, it will 
converge to stationary distribution vr. Here 7r(f) = d(v)/{2m), where d{v) is 
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the degree of vertex v. We often write Tr{v) as tTi,. The eigenvalues of P{G) 
are Aq = 1 > Ai > • • • > A„_i > —1. Let Amax = niax(Ai, |A„_i|). The rate 
of convergence of the walk is given by 



(6) |P^(x)-7r,| < (7r,.K)i/2A 



i 

max" 



For a proof of this, see, for example, Lovasz [20]. 

To ensure that the walk is both ergodic and that Amax = Ai we make the 
chain lazy; i.e., the walk only moves to a neighbour with probability 1/2. 
Otherwise it stays where it is. This shifts each eigenvalue up by 1, and so 
Ai > A„_i = 0, and (6) still holds. 

Now define the product graph H = H{G,k) = {Vh^Eh) to have vertex 
set Vh = and edge set Eu = E^ . The vertices v of H consist of /c-tuples 
V = {vi,V2, ■ ■ ■ , ■Ufe) of vertices Vi 6 Vq, i = 1, . . . , /c, with repeats allowed. 
Two vertices v,w are adjacent if {vi,wi), . . . , {vk,Wk) are edges of G. The 
graph H replaces k random walks (t) on G with current positions Vi and 
starting positions Ui by a single walk W^(t). 

If 5 C Vh, then T = T{S) is obtained from H by contracting 5* to a 
single vertex 7(5*). All edges, including loops and parallel edges are retained, 
producing a multigraph. Thus iir(7) = dniS), where dp denotes vertex 
degree in graph F. Moreover F and H have the same total degree (nr)^, and 
the degree of any vertex of F, except 7, is r^. 

Lemma 4. Let G be typical and let S C Vh be such that dniS) < 
j^2^k-i^k _ Lei F = G,H,r. Let be a lazy random walk starting at 

u £ Vp and let Tp be such that, for graph F = {Vp,Ep) and t > Tp, 



(7) 



I Dt^ N I ^ min^6y^(7ra;) 
max \P,,[x) — VTa; < . 



Then for k < n, 

Tg = 0(ln?i), Tff = O(lnn), and Ty = 0{k\Q.n). 

The above is a slightly modified version of the same lemma proved in [8] . 

We analyse our walks in the product graph and assume that we keep the 
chain lazy for the duration of the mixing time. At this point it is mixed, and 
we can stop being lazy. A lazy walk in the product graph maps to a process 
where all the walks move or don't move together. That is, with probability 
1/2, each walk independently takes a random step, and with probability 
1/2 none of them do. Consider the following two conditions: (i) interactions 
can only take place upon moving to a new vertex, and (ii) can only be 
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decreased (by 1) upon moving to a new vertex. It is not difficult to see (e.g., 
through coupling) that the laziness of the walk cannot affect the infection 
outcomes. Laziness affects time, but only during mixing periods since we do 
not keep the chain lazy thereafter. 

7.3. First visit lemma. Let G be typical, S C Vh be such that dniS) < 
f^2^k-i^k g^j^^ k < rf for sufficiently small e. Let F G {G, H, F} and let 
T = Tf. Let be a walk started at u G Vp. For t > T let ft = ft{u^v) be 
the probability that the first visit of iov ^Vp in the period [T, T+1, . . .] 
occurs at step t. 

For the graphs we consider, the following was shown to hold in [8]. 
Lemma 5. Let 
(8) A ^ 



KT 

for some sufficiently large constant K . Let v ^Vp and let 

where = 9j. + 0{k'^n^^^^^ + (/clnn)"^^*^^)) and r2+(l) denotes an arbi- 
trarily large constant. Then for all t >T, 

(10) ftiu^v) = (1 + 0{Tn,))-—^^ + O (rvr,e-^*/2 

For t > T let Ai,{t) be the event that Wu does not visit v during steps 
T,T + 1, . . . ,t. Then we have the following corollary. 

Corollary 6. 

(11) Pr(X(t)) = + ^^^^-^^ , + O (rV.e-^*/^) . 



7.4. Visit to a set of vertices: Contraction Lemma. The following lemma 
formalises the notion that we can deal with a first visit (after the mixing 
time) to a member of a set S of vertices of a graph H by contracting 5 into a 
single vertex 7 = 7(6") and instead deal with a first visit to 7 on this altered 
graph. 

Lemma 7 ([7]). Let W,,^ be a random walk in H starting at u ^ S, and 
let be a random walk in F starting at the same vertex ti 7^ 7. Let T be a 
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mixing time satisfying (7) in both H and T. Let Ayj{t) he the event that no 
visit was made to w at t >T. Then 



where the probabilities are those derived from the walk in the given graph. 

7.5. Which vertex in the set S was visited?. Given that we know a first 
visit has been made to a set S, we need to know the probabihty that the visit 
was made to a given v £ S. By applying Lemma 5, Corollary 6 and Lemma 
7 to the relevant vertices and contracted vertices, the following lemma gives 
the answer. 

Lemma 8. Under the assumptions of Lemma 5, let = Bi,{t) be the 
event that the first visit to S <Z H after T = O(fclnri) occurs at step t, and 
that this visit is to vertex v £ S. Let j{S) be the contraction of S, and let p^ 



be defined in T{S). Let 5 be the contraction of S\{v} and ps thus defined in 
T{S \ {v}). Let Py — PS = Pt>(l + e^). Let t > 2{T + L), where L = , then 



Proof. It is enough to prove the lemma for a two vertex set S = {u, v}, 
as one vertex can always be a contraction of a set. Let t be expressed as 
t = 2T + L + s, where s > L. Divide [0, t] into successive intervals of length 
r,s,r,L respectively, i.e. [0,T-1], [r,s+r-l], [s+T, s+2T-l],[s+2T,t]. 

Let A be the event that W(cr) {u, v} for cr G [T, s + T — 1] and that 
W{t) = u, but >V(fT) / n for cr G [s + 2r,t - 1]. Contract S and apply 
Lemma 5 to 7(5") in [T,T + s — 1]. Next, starting at s + T, apply Lemma 5 
to u for a visit at t. Thus 



Let Bu be the event that W{t) = u and W{cr) {u,v] for a £ [T,t - I]. 
Then fi„ C A and so Pr(fi„) < Pr(^). It follows from (13) that 




(12) 



Pr(E„) = (1 + 0(L7r^) + 0(e„))p,(l - p^f. 



(13) 



Pr(^) < (1 + 0(r7r^))(l - p^y (1 - p„)^p, 



(14) 



Vv{Bu)<Pu{1-P^)\1 + 0{Lt:^)). 



However, by contracting 5 we have that 



Vt{Bu U B,) = (1 + 0{Ttt^))p^{1 - p^)\ 
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and so 

Pr(i3„) > PriBuUB,,) -Pr (Bu) 

= (l-0(L7r^))(p^-p„)(l-P7)* 
(15) = (l-0(L7r^) + 0(e,))p,(l-p^)*. 

The result follows from (14) and (15). 

□ 

7.6. Particle pair meetings. Consider a pair of particles {x, y) (or (y, x)- 
we remind the reader that it is unordered, as per notation for graph edges). 
X and y being incident in G - that is, at the same vertex in G - maps in the 
product graph to a set of vertices S = {v = {vi,V2, ■ ■ ■ ,Vk) ■ Vx = Vy} C Vh- 
We can, therefore, calculate the probability for an xy meeting in G at time t 
by calculating the probability for the single random walk on H visiting 
S at time t. This in turn is done by the contraction described above, and 
calculating the probability of the walk on T visiting 7(5) at time t. 
These two times asymptotically equal by Lemma 7. 

More generally, let A C {(x, y) : x,y £ V,x y} he a set of particle pairs. 
If J^Ait) denotes the event 

{For some {x, y) £ A there is an xy meeting at time t}, 

in the product graph setting 7W^(t) maps to the event {W^(t) G S} where 
S = {v = {vi,V2, . . . ,Vk) ■■ Vx = Vy for some {x,y) £ A}. 

To use Lemma 5 and corollaries for the walk on T visiting 7(6*), we 
need to calculate the relevant (equation (9)). 

Lemma 9. Let A be a set of particle pairs. Let S{A) C Vh be such that 
V = {vi,V2, . . . Vk) £ S iff, for some pair (x, y) £ A, Vx = Vy where Vx (resp. 
Vy) is the position of x (resp. y) in G. Then inV = T{S{A)), 



\A\ f ^ /Plnn 1 



rn\ V (A;lnn)f^+(i) 

where 7 = 'y{S{A)), if k <rf and e is small enough. 

Proof. Let N = \S{A)\. A particular pair {x,y) £ A can be on n possible 
different vertices in G, and for each one, the other particles can be on n^~'^. 



imsart-aap ver. 2011/11/15 file: ViralProcessesl7.tex date: January 20, 2013 



VIRAL PROCESSES ON RANDOM REGULAR GRAPHS 13 

Thus < I A|n''~^. Further, N > N' where A^' is the number of G if such 
that only one of the particle pairs is incident. 



\A\'n!'-^ -0[\Ay'n 
I n 



\A\n^-'^ ( 1 - O ( 



The term (''^')n'^ ^ represents two pairs coincident at the same time, and 
(''^')n*^~^ represents more than 3 particles coincident at any vertex. Thus, 
N = lAjn'^"^ ^1 — O (^tt)) ^"^^ since each vertex of H has degree and 
contraction preserves degree. 



ri/^r'^ \ J J n \ \^ J J 



and 



|A|fcln?i\ / r-, f ^"^W ^{k^l'a 



Ttt^ = 0[ ^ 1 ( 1 - O ( — 1 1 = O 



n J \ \ n I I \ n 



since |^| < (2). 



Hence (9) becomes 



+ 0(A;2n-f^(i) + (A: lnn)-^+(i)) {l + O 



n V V n^(^) (A;lnn)f^+(i) 



□ 

That the ^^+(1) term can be an arbitrarily large constant is demonstrated 
in the derivation of Ry in [8] Lemma 19. 

During the mixing time, we cannot be precise in general. However, if the 
particles begin in general position with respect to each other then whp, the 
particles don't meet in the mixing time. This is seen in the proof of Lemma 
14. 

We say a visit to vertex u, or a particle-pair meeting is T-distinct if it 
occurs at least T steps after a previous T-distinct visit/meeting, or the start 
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of the walk. As suggested by the above, we can view the dynamics within 
the mixing time separately to the dynamics after the mixing time, up until 
the visit to a vertex or a particle-pair meeting. 

Lemma 10. Let A of particle pairs and let {x,y) G A (note, {x,y) is 
unordered). Let B(^^ y-^{t) he the event 

{no ab meeting for any (a, b) £ A in the period [T, t — 1] and xy meeting at t} 
Then 

(17, Pr(e,,,) = (l + o + j^^-L^) ) ± 

' \A\ f ( kHnn 1 \\y 

Proof. We use Lemma 8 with Lemma 9. Referring to equation (12), 
Lvr-y = 0(A;^(ln?7,)^/n) and 

_ / k^lnn k"^ \ _ f k''\nn 1 

(fclnn)f^+(i) J " V"^i^^ (/fclnn)f^+W 

where we have used the fact that the ^^"''(l) term is arbitrarily large in the 
RHS. Then 

Pr(%,„,) = (l + O + p-L^) ) 

1 ^fk^lnn ^ 1 



n V V n^^^^ (fclnn)^+W 

1^1 ^fk'^hin ^ 1 



OrU \ \ n^^W (A;lnn)^+(i), 

□ 

For notational convenience, we shall use p to stand for expressions of the 
form 

_ J_ / ( k^jhinf 1 

^ ' ^ ~ 6rn V V n^^^) ^ {klnn)^+W 

whence (17) can be more succinctly expressed as 
(19) Pr(B(,,,))=p(l-|A|p)* 

Observe that (1 - \A\pY = (1 - |A|p)*-i(l - 0{k^/n)), since \A\ = 0{k'^). 
This in turn has the form p{l — \A\pY~^ , which, when \ A\ = 1, resembles a 
Geom[-^-^) random variable taking value t. This last observation is exploited 
in the form of the A graph, described in section 8. 
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8. Case p = 1, ^ = oo. Two key ideas we use are that (1) the meeting 
times between pairs of particles are almost independent, and that (2) the 
meeting time for a pair {x,y) is roughly distributed as Geom{-^). It is 
self-evident by the independence of the random walks that xy meetings are 
independent of wz meetings. This is taken further by the assertion of (1) is 
that meetings between any of the (2) particle pairs are almost independent. 
Thus, consider the graph of the (2) particle pairs: if we weight each edge of 
this complete graph with a meeting time, then (1) and (2) imply the edge 
weights are (roughly) i.i.d Geom{-^-^). We use this intuition to define in 
section 8.1 a reduction of the viral process (for the SI case, i.e., ^ = 00), to 
a weighted graph on the particle set that captures the infection time of each 
particle. We formalise (1) and (2) in section 8.2. 

8.1. Interaction graph. The interaction graph T = (Vxi E-x) is a weighted 
complete graph on the particle set V, thus Vx = and E-^ = {(x, y) : x,y G 
V,x ^ y}. For a particle x, let t{x) be the time at which x is infected, 
or 00 if it never gets infected. For an interaction edge e = {x,y) G Er, 
let t(e) = m.m{t{x),t{y)}. Then the weight W"f{e) of the edge is a random 
variable defined as 
(20) 

wr{e) = min{t — t{e) : t > t{e) and there was an xy interaction at step t}. 

If the other particle was susceptible at the time of interaction, it becomes 
infected at the point of interaction. Therefore, T can be represented by a ran- 
dom (2)-vector {wr {ei) , wr {e2) , ■ ■ ■ ,wr{e^k-^)) where the edges are labelled 
in some arbitrary order. The distribution on T is the joint distribution of the 
edge weight random variables wrisi). An outcome T = z = (zi, ^2, . . . -^^fcp 

is a specific set of values for each random variable wx{ei). 

For a pair of particles x, y and some weighted graph F on the particle 
set, let DF{x,y) represent the weighted distance between x and y in F. Let 
dF{x) = m.m{DF{s,x) : s G 2^(0)}. Furthermore, for an edge e = (x,y), 
let dpie) = mm{dF{x),dF{y)}. The interaction graph allows us to relate 
weighted distance to infection time: 

Lemma 11. For a particle x, t{x) = df{x), and so for an edge e G Ej;, 
t{e) = dr{e). 

Proof. Consider a particle x. If x G X(0) then t{x) = = Dr{x,x) 
and the claim holds. Assume, then, that x ^ 2^(0). For any s G X(0), let 
P = (Po = s,pi,p2 ■ ■ ■ ,pi = x) be a weighted shortest path from s to x 
in T. Thus, Dr{s,x) = D{P), the weighted distance of P. For < i < 
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I - 1, observe < t{pi) + wr{pi,Pi+i) = t{pi) + Dr{pi,Pi+i) (since 

P is a shortest weighted path). Thus, if t{pi) < Dr{s,pi), then < 
Dris,Pi) + DriPi,Pi+i) = Dr{s,pi+i). Since t{s) = Dris,s), we conclude 
t{x) < Dj;{s,x). Since s was an arbitrary member of X(0) we deduce, t{x) < 
dr{x). 

Now, for X to have been infected, there must have been a series of infec- 
tions from one particle to another starting from some particle a € 2^(0). Let 
Q = {a = qo,qi, . . . ,qi' = x) be this series of infections. For < j < I' — I, 
since gj+i is infected by qj, we have = t{qj) + wr{qj,qj+i). Thus, 

iterating, we have t{x) = Yl^jZo ''^Tiqj , Qj+i) > Dr{a,x). Thus, since a was 
some particle in X(0), we deduce that t{x) > df{x), the shortest distance 
over all particles in I{0). □ 

In subsequent sections, this reduction is used to determine a completion 
time for the process and to derive another graph that gives the number of 
infected particles in the case ^ < oo. 

8.2. Interaction Graph Approximation. In this section we calculate the 
probability of a particular interaction graph. We claim that the joint dis- 
tribution of edge weights, is well approximated by the joint distribution of 
random edge weights of another complete graph on "P, A = {V\,E/^). The 
weight w\{e) of each edge e G E\ is an independent random variable with 
distribution Geom{q), where q = Recall = 1 when p = 1- 

Labelling the edges in A as with 1 < i < (2) , suppose z = {zi, Z2, ■ ■ ■ z/k\). 



Observe that with probability one, neither graph has an edge of infinite 
weight. Assume the edges of T and A are labelled in the same order. We 
shall prove the following 

Lemma 12. 



We say an edge e is active in the period [t{e),t{e) + wr{e) — 1]. The 
event {T = z} can be restated as follows: {For each = {x,y), the first xy 
interaction after time t(ej) is at time t{ei) + Zj}. 



Then 



(21) Pr(A = z) = Pr f\ WA{ei) 




Zi-l 




(22) 



Pr(T = z) = (1 + o(l))Pr(A = z). 
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By Lemma 11, this can be restated as 
(23) 

{T = z} = {For each edge = there is no xy interaction in the period 

[dr{ei) + l,dr{ei) + — 1] and there is an xy interaction at time df{ei) + Zi}. 

Let Aq be the set of active edges at time tq = 0. For example, in the 
case X(0) = {s}, Aq = {{s,x) : x £ V,x ^ s} whence |^o| = k — 1. Let Tj, 
1 < i < R for some R, denote the times at which the set of active edges 
changes. Let Ai be the set of edges active in the period [rj,rj+i — 1]. We 
call the Tj's epochs. A particular edge = {x,y) will become active at some 
epoch r = dx{ei) and remain so until r' — 1 (inclusive), where r' = dr{ei)+Zi 
is some epoch after r. Its period of activation may, therefore, cross a number 
of epochs, in which case it will be a member of a number of edge sets Aj. 

We analyse each period [tj + l,Tj+i] individually through application of 
(19). In the product graph H, we identify the set of vertices S{Aj) C Vh 
such that V £ S{Aj) iff for some (x,y) G Aj v represents x and y being at 
the same vertex in G. Thus, ii v = {vi,V2, ■ ■ ■ ,Vk), then Vr = Vg where r 
and s are the vector indices for x and y respectively. Since each pair is not 
to interact in the period [tj + l,Tj+i — 1], and since we assume p = 1, the 
walk W on must avoid the set S{Aj) until Tj+i — 1 (inclusive). Then, at 
Tj+i it must visit Sj C S{Aj). v G Sj iff it represents any of the {x,y) E Aj 
that interact at time Tj+i being at the same vertex in G. We deal with (non- 
)visits to S{Aj) and Sj by applying (19) to their contractions 'y{S{Aj)) and 

W do not prove (22) for all possible z, only those that are termed good. It 
will be shown that T is good whp. To prove this, we require the following 
minor adaptation of Lemma 20 in [8] . 

Lemma 13. Suppose the particles start walking on G (typical) with min- 
imum separation 

d = Q(lnlnn + Ink), 
and suppose r = 0{T^), then if e is small enough, 

(24) Pr(a given pair of particles x,y meet during r) = 0{t'^ /{r — V)'^^^). 

Let T = 0{k\n.n) be a mixing time and let L = T^, and let i = 2(T + L). 
The following is more than is required for the proofs of this section, but will 
be needed for the case p < 1. 
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Lemma 14. With high probability, T is good, which we define as a graph 
satisfying the following: (a) None of the (2) interactions that form the edges 

of T take place within l + T steps of each other, (h) (i) z-i < k'^nlnii 

and (ii) Zi < n^'^ for each i. 

Proof, (a) By an occupied vertex, we mean a vertex visited by at least 
one particle at that time step. We consider what happens during the first 
mixing time, when the particles start from general position, and also the 
separation of the occupied vertices when a meeting occurs. 

Suppose two (or more) particles meet at time t > T = 0{klnn). The 
position of the (at most) k particles under consideration, including the coin- 
cident pair, is close to the stationary distribution. There are 0{r'^) vertices 
within distance d of a given vertex of G, so the probability pa some pair are 
within distance d at any coincidence is 0{k'^r'^ /n). 

Now letting d = a(lnlnn + In A;) for some constant a, when a particle 
pair X and y have an interaction for the first time after the first of them 
is infected, then with probability pd = 0{k'^r'^ jn) there are some pair of 
particles within distance d of each other. Taking the union bound over all 
(2) pair-wise interactions, this is 0{k^r'^/n) = o(l) if ae is sufficiently small. 

By Lemma 13, taking the union bound over all pair-wise interactions 
again, the probability that there is a meeting during any of the (2) post 
meeting times i is 0{k^{i + Tf/{r - if/^) = 0{k^°lnn/{r - l)"^/^) = o(l) 
for sufficiently large a and sufficiently small ae. Hence this part of the lemma 
is true if e is small enough. 

(b) Consider a pair of particles x, y with x the first to be infected. Using 
Corollary 11, it is straightforward to show that E[/((a;— )-y)] = 0(n). For (i) 
sum over all (2) edges; the result follows by Markov's inequality. For (ii) 
use Markov's inequality and take the union bound over all the edges. The 
result follows by part (a) of the lemma. □ 

For part (a) of Lemma 14 we note that for a given a, ae can be made 
sufficiently small by making e small enough (as it is a constant with not 
restriction on how small it can be). We now prove Lemma 12. 

Proof of Lemma 12. We shall assume the z that is drawn is good. 
Recall that R is the number of times the set of active edges changes. Lemma 
14(a) implies R = (2) and so we only need to consider cases of z where, for 
each j, Sj is a single particle pair. It also implies that during the period 
[Tj -\- l,Tj -\- i], we do not need to account for interactions. By the Markov 
property, we can apply Lemma 10 sequentially to each period [tj + l,Tj+i] 
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and multiply the probabilities, giving (using the short-hand (19)) 

R 

(25) Pr(T = z) = JJp(l - \Aj\pY^+^~^^ 



j=0 



Let Aj = Tj+i — Tj. Then 



Pr(T = z) = p''llexp{-{l + 0{\A,\p))\Aj\pA,} 



j=0 



p^ exp < 



R 



-il + 0{k'p))pY,\Aj\A, 

j=0 



Furthermore, 
R 



R 



j=0 



j=0 eeEr 



and 



R 



R 



E ^{eeA,}Aj = E l{eeA,}(Tj+l " ^j) = ^e, 
j=0 j=0 

because this sums over all intervals [Tj,rj_(.i] in which edge e is active, that 

sum being Ze- Hence, J2f=o = 'Ee&Er = Ei=i 
Thus, 



Zi. 



Pr(T = z) = ^(a) exp < 



(^) 

-{l + 0{ep))pY,- 



i=l 



Now from (21), we have 
(^) 



Pr(A = z) = JJg(l - qY'-'^ = qi^) exp < 



(2) 

-(l + 0(l/n))g^, 



Therefore, 

Pr(T = z) 
Pr(A = z) 



P 



exp < 



j=i 



{{l + 0{l/n))q-{l + 0{k^p))p)Y,- 



1=1 
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- '1 + + ' ' 



nil) (A;lnn)^+(i) 



(26) = 1 + + 



n 



and 



(1 + o(i/„)), - (1 + o(.V))P = -o + 

By Lemma 14(b), whp 



n 



1 \ 

' n(fclnn)f^+W J ^''^ " n^W '(A;lnn)f^+W 



/A;5(lnn)3 fc7(lnn)^ k'^lnn 



Since the ^^"''(l) term in (26) and (27) is an arbitrarily large constant, 
they are 1 + o(l) and o(l) respectively, when k <rf for sufficiently small e. 
Thus, we conclude Pr(T = z)/Pr(A = z) = 1 + o(l). □ 

9. Case p = 1, ^ < oo. To analyse the SIR model, we build an inter- 
action graph but we modify the process by considering two phases. Phase 1 
assigns edge weights as before: an edge e = (x, y) is weighted by the time it 
takes for an xy interaction after one of the particles has become infected. It 
is possible, of course, that neither ever get infected when ^ is finite, in which 
case the edge (x, y) does not become active (and is not given a weight) in 
this phase. Phase 1 ends at Tend when there are no more active edges. At this 
point, there remains a set of susceptible particles S{Tend) that were never in- 
fected in phase 1. Phase 2 begins with an arbitrary particle in S{Tend) being 
given a pseudo-infection with an infinite infectious period, and we proceed 
with an SI process on the particle set S^Tend), giving them edge weights, 
as per the normal SI process. A pseudo-infection is the same process as 
an infection but we do not count pseudo-infected particles as infected; only 
particles infected in phase 1 are counted as infect. Phase 2 exists to maintain 
the same probability space as in the SI model (see remark below). 

Call the resulting (complete) graph ^. We remove an edge e in ^ if 
W\i,{e) > ^. We use the mapping : N^^) g^,^ where Qk is the space of all 
2(2) graphs on the particle set. maps a vector z representing a weighted 
complete graph to a graph G, where Cj € Eq iff zi < ^. 
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Lemma 15. In the SIR process a particle y is infected if and only if y G 
C{x) for some initial infective x G 2^(0), where C{x) denotes the connected 
component of x in f^{^)- 

Proof. Let G = f^{^)- y gets infected if and only if there is a path 
from y to some x G 2^(0) in G: If y got infected, there must have been 
a series of infections {x = pQ,pi, . . . ,pi = y) for some x G 2^(0). In each 
in this sequence we must have < ^, therefore the path 

{x=pQ,pi,...,pi = y) exists in G. 

Conversely, consider a path {x = pQ,pi, . . . ,pi = y) in G. Suppose pi 
was infected. Then either pi+i was infected before pi or pi was infected 
first. In the latter case, pi must have passed the infection to Pi+i since 
{pi,Pi+i) G G =^ w^{pi,pi+i) < ^. Since x = po was infected, inductively 
the whole path was. □ 

Remark We use this two-phase approach because wish to maintain the 
same probability space on ^ as in the SI model. Doing so makes it straight- 
forward to model the structure of f^{^) as an Erdos-Renyi random graph by 
removing edges with have weight greater than One may wonder why we 
do not employ the more obvious strategy of letting the infections run as an 
SI process then removing edges with weight greater than ^. Unfortunately, 
this approach fails, as can be demonstrated by the example in the appendix. 

We show that Lemma 12 holds for ^. 

Lemma 16. For z good, 
(28) Pr(^ = z) = (l + o(l))Pr(A = z). 

Proof. Our approach follows that for the SI model, except now we have 
to consider the two phases in turn. We are given particular graph weighting 
z that is assumed to be good and need to determine Pr(^' = z). Consider 
phase 1: Lemma 11 holds for all infected vertices and, as per the SI model, 
we have a sequence {Tj,Aj) defined by the edge weights. The active edge 
sets Aj are as before: an edge e = (x, y) is in a set of active edges Aj for the 
first time tj iff the first time x ox y became infected was at step Tj. Some 
particle pairs never become infected in phase 1; these will not be placed into 
any active set Aj in phase 1. Phase 2 is an SI phase but we only consider 
the set of susceptible particles after phase 1. 

The application of Lemma 10 is by the same means as for Lemma 12, 
giving probabilities Pr(^ = z; phase 1) and Pr(^' = z; phase 2) for phase 1 
and phase 2 respectively. For each phase, the proof of Lemma 12 applies in 
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a similar way as it does for the SI model giving 

Pr(^' = z; phase 1) = (1 + o(l))Pr(A = z; phase 1), 
Pr(^' = z; phase 2) = (1 + o(l))Pr(A = z; phase 2). 

Here Pr(A = z; phase X) is the probability that A has the edge weights 
determined in phase X. By the strong Markov property, 

Pr(^' = z) = Pr(^ = z; phase 1) • Pr(^ = z; phase 2), 

and so 

Pr(^ = z) = (1 + o(l))Pr(A = z; phase 1) • Pr(A = z; phase 2). 

Since the active edge sets generated in phase 1 and phase 2 are mutually 
exclusive, 

Pr(A = z; phase 1) • Pr(A = z; phase 2) = Pr(A = z) 

and the Lemma follows. □ 

Let Gk,q, the space of Erdos-Renyi random graphs with edge probability 
q = 1 — {I — q)^, where q = for the special case p = xf: = 1. With slight 
abuse of notation, we show the following: 

Lemma 17. An event 21 that occurs whp in Q^^^ occurs whp in f^{^). 

Proof. Lemma 14 holds for ^, so 

Pr(^' = z) = 1 - o(l) 

z good 

and by Lemma 16 

J2 Pr(M/=z) = (l + o(l)) Pr(A = z). 

z good z good 

Therefore 

Y Pr(A = z) = l-o(l). 

z good 



imsart-aap ver. 2011/11/15 file: ViralProcessesl7.tex date: January 20, 2013 



VIRAL PROCESSES ON RANDOM REGULAR GRAPHS 



23 



Let be the indicator variable that is 1 if 2t holds for G and otherwise. 
For a graph G on k vertices, Pr(^fc ,j = G) = Pr(/^(A) = G). Therefore, 

= j;Pr(A = z)2l|^^(,„ 

z 

= o(l)+ ^ Pr(A = z)%^(,)| 

z good 

= o(l) + (l-o(l)) Pr(M/ = z)2t|;^(,)| 

z good 

= o(l) + (l-o(l)) j;Pr(^' = z)2t|;^(,)}. 

z 

Hence, if E[2l|g^ = Pr(^ = z)2t|;^(,)} = E[2t|;^(^)}] = 

l-o(l). □ 

Proof of Theorem 1 when p = -(/; = l. For an Erdos-Renyi random 
graph space Qn,p on n vertices and edge probability p, the connectivity results 
are well known (see, e.g., [17]). By Lemma 17, whp statements on sizes of 
connected components in Qk^q carry through to f^{^)- Since Theorem 1 deals 
with the case of one initial infective s, the infected particles will be those 
in C{s) by Lemma 15. Cases (i) and (iii) of Theorem 1 are straightforward 
to see from the above. For case (ii), there is a unique giant component g of 
size (1 + o{l))Gk in f^{^). By the symmetry of the model, the probability 
any particular particle x being in g is \g\/k = (1 + o(l))C Thus, the giant 
component is infected with this probability. Otherwise, s will be placed in a 
component of size at most 0(lnA:). □ 

10. Two-particle System. We now address the case p < 1. To do so, 

we introduce the techniques we shall use by giving an informal treatment 
for a system with only two particles. In subsequent sections, we will give 
a more rigorous analysis for k particles, where, as per the assumptions in 
section 4, k < n^, where e is a sufficiently small constant. 

Let X and y be the two particles, with x being the initial infective and 
y being susceptible. Suppose that the assumptions stated above hold. We 
allow ^ < oo and/or p < 1. The former conditions means that y may never 
get infected, the latter condition means that it may take more than one 
meeting between x and y before an interaction takes place. Note, that if x 
and y were at the same vertex at time t, and happen to move to the same 
neighbouring vertex in the next step, then this counts as another meeting. 
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By (19) with A = {(x,y)}, 

(29) Pr(^(,,,))=p(l-p)*-\ 

Now, suppose X and y have just stepped to the same vertex v. With prob- 
abihty p there wih be an interaction. After this, they will move again, either 
to the same neighbour of v with probability 1/r or to different neighbours 
with probability (r — l)/r. Let 

4> = Pr(No xy interaction before they move apart) 



S"-"(r(-;) 



1\ (l-p)(r-l) 



1 + P 



The following lemma is from [8]: 



Lemma 18 ([8]). Let G be a typical r -regular graph, and let v be a vertex 
of G, treelike to depth Li = [ei log^ nj . Suppose that at time zero, two inde- 
pendent random walks (Wi,yV2) start from v. Let {a,b) denote the position 
of the particles at any step. Let S = {{u,u) : u £ V}. Let f be the probability 
of a first return to S within T = 0{klnn) steps given that the walks leave v 
by different edges at time zero. Then 

Using this lemma, let 
(pT = Pr(No xy interaction before being apart more than T time steps) 

(30) ~ E^^/^"'(i-/) = T^- 

Now, assuming x and y start at the same vertex, 

Pr{xy interaction occurs within T time steps) 
~ Fr(xy interaction occurs before x and y have been apart more than T steps) 

(31) =i-,, = i-i(i_Z).i^^=,. 

Recall ip was defined in (1), and observe that p < ijj < 1 with ^ = 1 iff 
p = l. 
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Clearly, the number of T-distinct meetings in [0, t] is at most t/T. Subject 
to slightly more stringent restrictions, it can be shown that we have 

(32) Pr (there are i T-distinct meetings in [0,t]) ~ ^.^p*(l — 

i.e., it is approximately distributed Bin{t, g^)- The probability that there 
are no interactions in any of the i intervals [tj , tj + T] where tj is the time 
of the j'th T-distinct meeting is (1 — ipy. Thus 



/A 

Pr(there are no interactions in the period [0,t]) ~ ( . J {p{l — ip)y{l — pY * 

i=0 

~ (1-Vp)* 

Hence, 

(33) Pr(j/ gets infected within time ^) ~ 1 — (1 — ipp)^ . 

When p = ip = 1, (33) looks similar to the bracketed terms in (2). This is, 
of course, not a coincidence since the bracketed term in (2) is essentially the 
probability that an infection is passed between a pair of particles if one of 
them had been infected, and therefore, cp is effectively the expected number 
of other particles that are infected by a particular particle. 

11. The case p < 1. In this section, we formalise some of the ideas of 
section 10, extended to k particles. 

As per section 8.1, we deal with sets Aj. However, in this context, we say 
there is no interaction for any pair {x,y) £ Aj in the period [tj , tj^i — 1] then 
a particular pair {x, y) interact at time Tj+i. We shall justify the assumption 
of only a single pair interacting by extending the justification for the p = 1 
case. We work in the product graph framework in similar fashion as for the 
p = 1 case. 

As was the case for the two-particle system, p < 1 implies there may have 
been a number of meetings before an interaction took place. Of course, a 
meeting may have occurred for any pair in Aj, not necessarily only the pair 
that interact at time Tj+i. However, when some pair do meet, they do not 
interact. 

For notational convenience, we shall write the period [rj,Tj+i] as [0,t]. 
Within this period and before the final interaction, there will be some num- 
ber of meetings in general. In the first T steps after a (non-interactive) 
meeting, there is no interaction with probability 1 — ijj' , where ip' is given in 
Lemma 19. We then continue the walk for 2(T + L) steps unobserved, after 
which it is observed. We redefine i, to he i = 2{T + L) + T = 3T + 2T^. 
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Lemma 19. Suppose G is typical and that at time zero, two particles x, y 
positioned on a tree-like vertex v interact with probability p. Let ip' be the 
probability of an xy interaction in [0,T — 1]. Then 



n 



where 



r-2+p- 



Proof, x and y start at the same vertex at time t = 0, and interact 
with probability p at time 0. Suppose the first xy interaction occurs at time 
r > 0. Let 21 = {r < T}. Now suppose x and y are incident (i.e., at the 
same vertex) at times = to,ti,. . .. Suppose tM is the smallest tr in this 
sequence such that tM+i — tM > T. Hence, tM + T is the first time that x 
and y have been apart T steps. Let ^ = {t < tM + T}. 

We first demonstrate that Pr(2t) ~ Pr(*B). 

Pr(^) < Pr(2t) 

< PrC^) + Pr({T > T} n {r < tM + T}) 

< Pr(^) + Pr({T > T} n {x, y meet in [T, 2T]}) 

< Pr(^) + Pr({x, y meet in [T, 2T]}) 



Using (11), 

Pr({x, y meet in [T, 2T1}) = 1 (1 + 0{Ttt^)) _ ^ - 

{l + {l + 0{T7r^))7r^/R^f^ V ^ 

= 0{{klnnf/n) 

Thus 

Pr(^) < Pr(2l) < Pr(^) + 0{{klnnf/n), 

i.e., Pr(*B) = Pr(2l) + 0{{klnn)^/n) 

To get the correcting factor in (30), observe that the sum "^^yi 4>^ (1 — /) 
assumes every vertex at which the particles part is tree- like. Let W(t) be the 
walks on G and let X{t) be a walk on an infinite r-regular tree T rooted at 
the start vertex which is assumed to be tree-like. We couple W and X until 
time Li. Since G and T have the same structure out to Li, the two processes 
are identical until t = Li. Let Yj = dist(x, y) in G. It is shown in [8], proof of 
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Lemma 17, that for Pr(yi J < Li/2 = 0{n~^^^'^'^) where the IS an ar- 

bitrarily large constant. It was also shown that, subject to /c < n*^ for a suffi- 
ciently smaU e, Pr(the walks meet in [Li,T] and Y^i > Li/2) = 0(n~^(^)). 
Thus, 



1 , „ 

Now we include the error term for the asymptotic equality in line (31). 

(34) i_i(W) ^ P{r-l) (^^pf-m) 

^ ' l-4>f r-2 + p\ V 

= v(l + o(n-^W)) =Pr(!B). 

Note in (34) we have used the assumption that p is a constant to absorb 
(functions of) it into the O term. Thus 

Pr(2l) = Pr(^) - O (^^) = ^ (l + O (n'^W)) - O (^^) • 

p is a constant means ■0 is a constant and so defining ?/^' = Pr(2l), 

,/ , / ^ / f/clnn)^' 



□ 

We require the following lemma 

Lemma 20. For each edge e £ Zx, wy{e) > ^Jn whp. 

Proof. We use Corollary 6. For a set of particle pairs A, let S = S{A) 
and 7 = 7(6*). Recall that A'y{y/n) is the event that the random walk in the 
contracted graph T does not visit 7 in the period [T, T + 1, . . . , y/n\. Recall 
also that = Or + o(l). Hence, 

Pr{A,iV^)) = + ^^^"^^^ ^ + O (T%e-'^/') . 

{l + {l + 0{T7T^))7r^/R^)^ ^ ^ 

= (1 + 0(T^^))(1 - 0{^/^7T^/R^)) + O (T%e-^^^/^^ 

= l-Oik^n). 

Therefore, the probability that any of the (2) edges has weight less than ^/n 
is 0{k'^ /n) = 0(1) for small enough e. □ 
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By Lemma 20 we have t > ^/n > T . For the period [0, t], 

t 

Pr(first interaction at t) = ^^Pr(z meetings and first interaction at t) 

i=0 

= Pr(i meetings and first interaction at t) 

i=0 

t 

+ Pr(i meetings and first interaction at t) 

i=T2+l 

We can assume T > Inn, and so, 

t t 

Pr(i meetings and first interaction at t) < (1 — ipY 

< (1 - 

Thus taken over all (2) interactions, we still have 0(l/n^(''^")) probabil- 
ity of any interaction involving more than meetings. We have used the 
assumption that p is constant. 

In our analysis, we use the first visit lemma on the product graph for 
sequences of visits thats are ^-distinct, where a visit is ^-distinct if it is at 
least £ time steps after the last ^-distinct visit. Thus, in a period [0, i], the 
number of ^-distinct visits is at most t/i. 

We extend Lemma 14 with the following: 

Lemma 21. Consider all the meetings between active particle pairs over 
the whole process. If k < and e is small enough, then whp; (a) Only one 
pair meet at a time and no other pair meet within £ = 3T + 2T^ steps. That 
is, particle pair meeting times are i-distinct. (h) Any active pair that meets 
at some step r and does not interact in the period [r, r + T — 1] does not 
meet in the period [r + T, r + . 

Proof, (a) We use the same principles as in Lemma 14(a). In that proof 
we consider precisely (2) meetings; when p < 1 we let the random variable 
X > (2) count the total number of meetings. Note E[X] = ©((g)/^) = 
0{k'^) since ifj is constant. 

Recalling from the proof of Lemma 14(a) pd = 0{k'^r'^ /n), union bounded 
over all meetings this is p^X with expectation 0{k'^r'^/n) = o(l) for small 
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enough ae. Thus, when an active particle pair meet, each of the other par- 
ticles are distance at least d from each other. The second part of the proof 
as per Lemma 14(a) is adapted in the same way. 

(b) We can assume r = 0. Suppose a pair of particles x, y meet. Let 21 
be the event {xy interaction in the period [0, T — 1]}. Let £ be the event 
{x,y meet in the period [T, £]}. Then 

^ ' ^ ~ Pr(2l) 1 - ij' 

Thus, since il^' is a constant, Pr(£|2l) is of the same order as Pr(£). The 
rest of the proof is similar to part (a). □ 

Let T = (ti, tj, tj_|_i) be a sequence such that (i)X]sii*s = *i (ii) > 
e = 3T + 2T^ for each s G {1,2, . . . (iii) U + U+i > £ + T and (iv) 

< U+i < T. 

We use T to represent having i ^-distinct meetings in the period [0, t] 
with an interaction at step t. Specifically, we interpret r as representing £- 
distinct meetings at times ti, {ti + t2), {ti +^2 + ^3), • • • 1 Yll=i and the first 
interaction occurring at time t. Note that it may be the case that tj+i = 0, 
in which case X]s=i = t, meaning that the first interaction occurs at the 
i'th (i.e., final) ^-distinct meeting. 

We let the random variable H represent a history in [0, t], that is, a run of 
the process in the period [0,t]. The event {H = r} means that the history 
that occurs is represented by the sequence r as interpreted above. 

Observe tjj' = Pr(2l) = Yli=o '^^^sre pi = Pr (interaction occurs at step i). 
Hence po = p. For a given r, let pr = Pu+i- 

We can now show the following lemma. From this, the proof of Theorem 

1 is in the same manner as for the p = 1 case in section 9. 

Lemma 22. Pr (T = z) = (1 + o(l))Pr (A = z). 

Proof. We use Lemmas 5, 7, 9 and 10. Let pj = \Aj\p where p is from 
(18). 

Pr(H = r) =pril + 0{L7T,) + 0{ej))p{l-pjY^ 

(35) X n (1 - ((1 + ^(^^^■)) (TrfyiTT + 0{Tn,e-''^/')^ . 
We elaborate on the structure of this expression. The term 

(36) ((1 + 0(Tvr,))^^-^^^ + 0(rvr,e-^W2)^ . 
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is the probability of no pair in Aj meeting in — 1 time steps then a meeting 
of at least one pair. That is, if we define t* = X]h=i then the expression 
represents no meeting in the period [tl^i,tl — 1] followed by a meeting at 
time t* . 

Following the meeting at t*, there is a T period in which there is no 
interaction, thus giving the multiplication by the expression (1 — -0'). 

It should be noted that whilst (36) is the probability of a meeting between 
some (possibly more than one) pair of particles {x,y) in the set Aj, the 
expression (1 — ip') is the probability for a single pair not interacting within 
the mixing time T. Thus, we are implicitly assuming that only one pair 
{a,b) meet, which is proved to be the case whp in Lemma 21. Following 
that period of length T, we can ignore the unobserved i steps by Lemma 
21(b) for (a, b). By Lemma 21(a), the whole period T + £ can be ignored as 
far as other particles are concerned. 

The component 

(1 + 0{L7Tj) + 0{ej))p{l - p,Y^ 

is the probability of no meeting between any pair in the final inter-meeting 
period tj, then the particular pair (x,y) of Aj meeting at time t*. Since 
t* < t iS tj+i > 0, the probability of interaction happening at time t is given 

by Pu+i = PT- 

For a history t = (ii, ij, tj+i), let 

F(T) = F(ti,...,t,+i) = ^^l^^^. 

PU+l 



Fir) =(1 + OilTTj) + 0(6,))p(l - pjY 



^ (1 - ^I^T'pT' 

i-l 

X TT ( 1 + OiTTT,) +oi^] e-^*'/2+o(i.p.) 



n(i.o,r.,).o(^) 



Recalling that Ry is a positive constant we have = @{7ry). Assuming 
mins=i,...,i > r^, and using (8) then 



^ (1 - n-'pT' 

i + o(r7r,) + o(re-^'(^+«(^^»' 



F{T)={l + 0{L7r,) + 0{e,))p{l-pj) 

' I z — 1 ' 



•j-i 
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Also 

= (1 - pjf~' (1 + 0{tp]) - 0{ipj) + 0{Tp,)) 



Thus, 



F{t) ={l + 0{L7rj) + 0{ej)) 



X {i + 0{tp])-Oiipj) + 0{Tpj)) 
xil-ijT'pp^\l-p^Y~' 



We simplify the correcting factor above, pj = \Aj\p = 0('/rj), \Aj\ < {^) and 
ej = —O [k^ In n/n). We are also assuming i < T^. Hence 

1 + OiLTTj) + 0{ej) = l + 0{T^TTj) -0{k^lnn/n) 

= 1 + 0{k^{lnnf\Aj\/n) - O {k^lnn/n) 
= l + 0{k^{lnnf/n). 



1 + 0{tp]) - 0{ipj) + 0{Tpj) = 1 + 0{tk^/n^) - 0{T^k'^/n) + 0{Tk^/ 



n] 



1 + Oitk^n^) - 0{k\lnn)yn). 



i + o(rvr,) + o(re-^'(^+o(^^)) 



1 + 0{T\,) + O (T^e-^'('+0{. 
l + 0{k^(lnnf/n) 



Thus, we have 



Fir) 



1 + 



A;^(lnn) 



n 



+ 



tk^ 



p{{i-^l,')p^y-^[i-p^) 



\t-i 



Observe that for a given t, fixing ii, . . . and allowing tj+i to vary 
from to T, determines ti and gives a set of histories that we can sum 
together as follows. 
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5]Pr(H=(ti,...,t,_i,if\r)) 



r=0 



1 + 
1 + 



n 



n 



n 



+ ^ 



r=0 



where ^ = t - {h + . . . + ti_i + r). 

Note r satisfying conditions (i)-(iv) define all ^-distinct meetings in [0, t] 
except those for which ti + tj+i < i + T, that is, those for which t*_-^ G 
[t — {i + T),t — i]. This window of size T can be ignored since, from the 
stationary distribution, the processes has probability 0{k'^T/n) of falling 
into it. Taken over all (2) periods this is still o(l) for k < rf and e small 
enough. Thus, for an integer i > 1 let Hj denote the set of all r that satisfy 
conditions (i)-(iv) and have tj+i = 0. Then for the period [0, i], and assuming 
all meetings are ^-distinct as per Lemma 21, 

Pr (first interaction at t) 

= O f ^^^^^^ j + Pr(i ^-distinct meetings and first interaction at t) 



o 



o 



1 



j^f2(ln n) 



+ 



1 + 



1 + 
A:^(lnn)^ 



A;^(lnn) 



n 



+ 



tk^ 



^i,'p{{i-^p')p^y-\i-p^) 



\t-i 



n 



^ ^ ^ ^ i=l T£Si 
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Now 




where Si is a set of all sequences in which there is some tg < i for 1 < s < i— 1 
and for which ti < £ + T. 
Thus 



Pr (first interaction at t) 



1 + 



1 + 



k^{lnnf tk^ k^lnnY 



+ 



n 



k^{lnnf tk^ F(lnn) 



+ 



n 



4\ n 



(we assume all O terms are o(l)). 
Now, 

-^ = l + 0(p,) = l + of-') 
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i=l ^ ^ 



- E 

i=r2+i 

(1 - ip'pjY - ©(l/n^^'""")) 



Therefore 
(37) 

Pr (first interaction at t) 



1 + 



+ 
n 



^V(l - i^'Pj 



I \t-i 



Thus, assuming the O term is sufficiently small, we repeatedly apply (37) 
to calculate the probability of an interaction graph, as we did for the p = 1 
case. Recalling that y/n <t< n^'^, by Lemma 14 and Lemma 21, 



Pr(T = z) 



n 

j=0 



l + O 



k^{lnn)^ tk^ /c'^(lnn) 



4\ n 



+ 



n 



ilj'pil-ij'lAjlp)^^ 



= (l + 0(n-"">+tV(fclnn)"+W)) (^j JI{1 - ^'\Aj\p) 

The term indicates an arbitrarily large constant. 

Continuing in the same fashion as for the p = 1 case, 

(2) (2) 

f[il-^P'\A,\pf^ = nexp{-(l + 0(|^,|p))V''|^,-|pA,-} 

j=0 j=0 

e^p\-{l + 0{k^pWpY,Zi 

i=l 
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Letting q=-ff^, 



1=1 



Pr (A = z) = JJ q{l - qY' ^ = q(^) exp < 
Thus 

Pr (T = z) 



(2) 

-(l + 0(l/n))g^. 



i=l 



Pr (A = z) 



1 + 



1 1 



exp|-(l + 0(A:Vn))^VE£U.} 
exp|-(l + 0(l/n))gE£U4 



= exp ^ O I + Tpf— TT^+Tiy ) 7— X] 



exp<'-(l + 0(l/n))gE-!i^. 



^(1) (A;lnn)f^+(i) / OrJi ^ 
^ ' ' 1=1 



1 + 0(1). 



□ 



12. Completion Time. In [8] an expectation of ^^^In/c was deter- 
mined for the completion time of a broadcasting model on k particles that is 
equivalent to the SI model with /q = 1 and p = 1- Generalising to '^0- In /c 
for p < 1 is straightforward. To get a convergence in probability to the 
same value, we apply a theorem from [16]: Assign each edge («, j) of a com- 
plete graph on n vertices a random weight Yij. The weights are assumed 
to be independent and identically distributed, non-negative and satisfying 
Pr(lij < t) = t + o(i) as t — )• 0. Let Xij be the minimal total weight of a 
path between a given pair of vertices 

Theorem 23 ([16]). For any fixed i, 

(38) ra^^j<n X,^ 2. 

Inn/n 

We apply this below. 
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Proof of Theorem 2. The 6x^(1) distribution is valid for Theorem 23. 
Since for a random variable Y distributed as exp{X) we have Pr(Y < t) = 
1 - e"^*, then Pr{Y < t/X) = 1 - e"* = Pr([/ < t) where U is distributed 
as exp{l). Thus, we can multiply the RHS of (38) by 1/A for edge weights 
distributed exp{\). Thus, when A = ip/{9rn), then applying Theorem 23 to 
the complete graph on k vertices, we have, for any fixed i, 

n\nk/k 'ip 

We assume i to be the initial infective. However, an edge in A has 

Pr{wA{e) >t)=(l--^) = exp + 0(l/n))/^ 
\ Urn J I Urn 

Thus, if Y is distributed as exp{ip/{9rn)), then 
Pr(y>t) _ exp{-|^} 



Pr(u;A(e) > t) ^ o{l/n))^^j 



if t = o{n'^), which is the case for good Zx- Thus, 'PY{w\{e) < t) = 
(1 + o(l))Pr(y < t). Hence, if i is the initial infective, then whp, the last 
infection occurs within time (1 it e)^^^ In A: where e is any constant. □ 
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13. Appendix. We demonstrate a failing example as per the remark 
in section 9. Suppose there are four particles a, b, c, d with a being the ini- 
tial infective. Suppose the following {particlepair, timestep) meetings take 
place: (o6,9), (ad, 11), (6c, 18), (cd,22), (cd,27), (ac, 100), (M, 100). The 
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SI interaction graph would have the following edge weights: wriajb) = 9, 
wr{a,c) = 100, wria,d) = 11, wr{b,c) = 9, wr{b,d) = 91, wric,d) = 11. 
Therefore, if ^ = 10 and we had a rule to remove any edge with weight 
greater than this, it would leave d isolated, suggesting it does not get in- 
fected. However, it clearly gets infected by the chain a b c d. 
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